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Abstract 

Wc have searched for the decay loI^v (l = e ov fi) in 78 fb^^ of T{4S) data (85 mihion 

BB events) accumulated with the Behe detector. The final state is fully reconstructed using the 
CO decay into 7r+7r~7r°, combined with detector hermeticity to estimate the neutrino momentum. 
A signal of 414 it 125 events is found in the data, corresponding to a branching fraction of (1.3 it 
0.4 ± 0.2 ± 0.3) X 10"^, where the first two errors arc statistical and systematic, respectively. The 
third error reflects the estimated form-factor uncertainty. 

PACS numbers: 13.20.He, 14.40.Nd, 12.15.Hh 



3 



The magnitude of Vub plays an important role in probing the unitarity of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix The cleanest way to constrain this quantity is either 
by measuring the decay B XJi^ inclusively, or by reconstructing one of its exclusive 
sub-modes. As to the latter, the decay modes B ixlv and B — > plu have already been 
observed 0, 0|- In this letter, we present a study of the decay B^ ujI^p which has 
not been measured so far E . Bl • Using three different form-factor calculations, ISGW2 0, 
UKQCD ^ and LCSR jlO|, we extrapolate the decay rates to the full range of lepton 
momentum and measure the branching fraction of this decay. 

The analysis is based on the data recorded with the Belle detector |ll| at the asymmetric 
e^e~ collider KEKB [l^] operating at the center-of-mass (cm.) energy of the T(45') reso- 
nance. KEKB consists of a low energy ring (LER) of 3.5 GeV positrons and a high energy 
ring (HER) of 8 GeV electrons. The T(45') dataset used for this study corresponds to an 
integrated luminosity of 78.1 fb~^ and contains (85.0 ± 0.5) x 10^ BB events. In addition, 
8.8 fb~^ of data taken at 60 MeV below the resonance are used to study the continuum 
{non-BB) background. 

The Belle detector is a large-solid-angle magnetic spectrometer consisting of a three-layer 
silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of- flight scintillation 
counters (TOP), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) 
located inside a super- conducting solenoid coil that provides a 1.5 T magnetic field. The 
responses of the ECL, CDC (dE/dx) and ACC detectors are combined to provide clean 
electron identification. Muons are identified in the instrumented iron flux-return (KLM) 
located outside of the coil. Charged hadron identification relies on the information from the 
CDC, ACC and TOP sub-detectors. 

Pull detector simulation based on GEANT is applied to Monte Carlo simulated events. 
This analysis uses background Monte Carlo samples equivalent to about three times the 



integrated luminosity. The decay B D*!^ is simulated using a HQET-based model |14 
The ISGW2 model is used for the decays B Dlv and B D**lv. The modes B D^*^ttIv 
are simulated according to the Goity- Roberts model |15|]. The ISGW2 and the De Pazio- 
Neubert model |16i] are used to model the cross-feed from other decays B — »• XJ,v. 

Events passing the hadronic selection 17| are required to contain a single lepton (electron 
or muon) with a cm. momentum p* [l8| between 1.8 and 2.7 GeV/c. In this momentum 
range, electrons (muons) are selected with an efficiency of 92% (89%) and the probability 
to misidentify a pion as an electron (a muon) is 0.25% (1.4%) [Til li^. 

The missing four-momentum is calculated, 

Pmiss = PHER + PLER -^Vi , 

i 

Emiss = -EhER + EiER ~^^Ei , (1) 

i 

where the sums run over all reconstructed charged tracks (assumed to have the pion mass) 
and photons, and the labels HER and LER refer to the two colliding beams. To reject 
events in which the missing momentum misrepresents the neutrino momentum, the following 
requirements are apphed. The total event charge must be close to neutral: IQtotI < 3e; the 
polar angle of the missing momentum (with respect to the beam direction) is required to 
lie within the ECL acceptance: 17° < ^^miss < 150°; and the missing mass squared, m^j^g = 
-^miss~Pmiss' required to be zero within about ±3 standard deviations: |^«missl < 3 GeV^/c^. 
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For generic B Xulv events, the efficiency after applying these requirements is 11%, 
and the resolution in the magnitude of the missing momentum is around 140 MeV/c. As 
the missing energy resolution is worse than the missing momentum resolution, the neutrino 
four- momentum is taken to be (|pmiss|,Pmiss)- 

Pairs of photons satisfying > 30 MeV, p^o > 200 MeV/c and 120 MeV/c^ < 772(77) < 
150 MeV/c^ are combined to form 7r° candidates. The decay u 7r"'"7r~7r° is reconstructed 
using all possible combinations of one 7r° with two oppositely charged tracks. Combinations 
with a charged track identified as a kaon are rejected, and the following requirements are im- 
posed: p* > 300 MeV/c, 703 MeV/c^ < m(7r+7r-7r°) < 863 MeV/c^. The Dalitz amplitude, 
A oc |p7r+ X Pn-\, is required to be larger than half of its maximum value. 

The lepton in the event is combined with the uj candidate and the neutrino. To reject 
combinations inconsistent with signal decay kinematics, the requirement \ cos6by\ < 1-1 is 
imposed, where 

2E*E* - ml - ml 
cos Oby = \ , f , 2 



and E^, p*^ and are fixed to -Ebeam = V ErerEler, \I E*^ - m% and 5.279 GeV/c^, 
respectively. The variables Ey, Py and my are the measured cm. energy, momentum and 
mass of the Y = u + I system, respectively. For well-reconstructed signal events, cos^^y is 
the cosine of the angle between the B and the Y system and lies between —1 and +1 while 
for background the majority of events are outside this interval. 

For each B^ uol^v candidate, the beam-energy constrained mass Mbc and /S.E are 
calculated, 



^bc = ^{El^.rr.Y-\pl+pl+plW 

AE = {E: + E! + El) - El,,^ , (3) 

and candidates in the range Mbc > 5.23 GeV/c^ and |A£^| < 1.08 GeV are selected. On 
average, 2.5 combinations per event satisfy all selection criteria, and we choose the one with 
the largest uj momentum in the cm. frame. Monte Carlo simulation indicates that this 
choice is correct in 77% of the signal cases. 

In BB events, the two B mesons are produced nearly at rest, and their decay products are 
uniformly distributed over the solid angle in the cm. frame. Conversely, continuum events 
have a jet-like topology. We exploit this property to suppress continuum background with 
the following quantities (defined in the cm. frame): the ratio R2 of the second to the zeroth 



Fox- Wolfram moment [2l| which tends to be close to zero (unity) for spherical (jet-like) 
events; cos6'thrust, where ^thrust is the angle between the thrust axis of the ul system and the 
thrust axis of the rest of the event; and a Fisher discriminant |22] that selects events with 
a uniform energy distribution around the lepton direction. The input variables to the latter 
are the charged and neutral energy in nine cones of equal solid angle around the lepton 
momentum axis. The selection £s/(£s + ^^b) > 0.9 is applied, where £3 (£b) is the product 
of the signal (background) p.d.f.'s of these three quantities. This selection is 56% efficient 
for signal decays and eliminates 92% of the remaining continuum background. 

The signal yield is determined by a three-dimensional binned maximum likelihood fit 
taking into account finite Monte Carlo statistics j23^. We use nine 240 MeV wide bins in 
AE, eight 20 MeV/c^ wide bins of m(7r"'"7r~7r°), and three 300 MeV/c wide Pi bins. The signal 
resolutions in AE and m{n~^n~7i'^) are about 140 MeV and 11 MeV/c^, respectively. The 
backgrounds from the remaining continuum events and from BB events in which the lepton is 
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TABLE I: The result of the fit assuming ISGW2 form-factors for ujI^v. The uncertainties 

quoted are statistical only. 







range 


(GeV/c): 






1.8 - 2.1 


2.1 - 2.4 


2.4 - 2.7 


1.8 - 2.7 


raw yield 


16,777 


4639 


326 


21,742 


continuum 


234 ± 45 


294 ± 50 


78 ±26 


606 ± 72 


other bkgrds. 


211 ±9 


216 ±9 


28 ±3 


455 ± 13 


subtr. yield 


16, 332 ± 46 


4129 ±51 


220 ± 26 


20, 681 ± 74 


B+ ujI+v 


101 ± 31 


193 ± 59 


89± 28 


383 ± 118 


B XJu 


339 ± 151 


466 ± 142 


125 ± 19 


930 ±312 


B XJu 


15, 755 ± 238 


3592 ± 64 





19, 348 ± 289 


sum 


16, 196 ± 284 


4251 ± 166 


215 ± 34 


20, 662 ± 442 



misidentified or does not originate directly from a B decay are subtracted from the raw yield 
bin-by-bin. The continuum background is estimated using the off-resonance data (scaled to 
the on- resonance luminosity). The fake and non-primary lepton backgrounds (which account 
for only about 2% of the raw yield) are determined from the simulation. The signal yield, the 
background from B Xylv decays and the background from B —>■ XJu decays are fitted 
after this subtraction. In the region defined by 763 MeV/c^ < m(7r'^7r~7r°) < 803 MeV/c^, 
I Ai^l < 360 MeV the signal purity is almost three times higher than the average in the whole 
fitting space; the broader fit ranges in AE, m^Ti^Ti^TT^) permit a reliable determination of 
the signal and background components. The distribution shapes of the three fit components 
are taken from Monte Carlo simulation. The contribution of each component is described 
by a single parameter. 

Tableland Fig.^show the result of the fit assuming ISGW2 form-factors for B^ uol^v. 
We find 383 ± 118 signal events with a statistical significance of 3.8 standard deviations. 
The latter is defined as y^— 2 ln(>CB / -^^s+b) , where £s+b (-^b) refers to the maximum of the 
likelihood function describing signal and background (background only). In addition to well- 
reconstructed signal decays, the signal component of the fit also includes candidates in which 
the lepton stems from a B^ ujl^v decay but the uj failed to be reconstructed properly. This 
sub-component is shown separately in Fig. ^ It scales with the actual signal and amounts 
to 36% of the signal component within the 763 MeV/c^ < m(7r+7r^7r°) < 803 MeV/c^, 
\AE\< 360 MeV region. 

The fit is repeated for the three form-factor models considered for B^ —>■ ujl^v. For each 
model, the signal yield N{B^ —>■ ul'^u) is determined and the branching fraction B{B~^ 
iul^u) is calculated according to the relation N{B'^ u!l~^u) = N{B'^) x B{B^ ujI^u) x 
B{uj — * TT+TT'Tr") X (ee±e^), where N{B~^) is the total number of charged B mesons in the data 
(assumed to be equal to the number of BB events), B{lj 7r"'"7r~7r°) = (89.1 ± 0.7)% ,2J| 
and eg (e^) is the model-dependent selection efficiency for cjez/ {u^v) signal candidates. 
The results are presented in Table in] Averaging the central values and the statistical 
uncertainties over the three models (giving equal weight to each), a branching fraction of 
(1.3 ± 0.4) X 10~^ is obtained. The spread around this average value amounts to 0.3 x 10~^, 
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ujl~^iy after requiring 763 MeV/c^ < m(7r+7r"7r°) < 803 MeV/c^ and |A^| < 360 MeV, respectively. 
The data points are the background subtracted yields. The open, hatched and doubly-hatched 
histograms correspond to the ujl~^v signal, the B — > X^lv background and the B — > 

Xclv background, respectively. In the highest momentum bin, the contribution of signal candidates 
in which the lepton stems from a B^ ujI^u decay but the uj is reconstructed improperly is shown 
by the dashed histogram. 



and provides an estimate of the form-factor model uncertainty. 

The experimental systematic error is 18.1% of the branching fraction (Table IIII|) . or 
0.2 X 10"'^ in absolute. The largest contribution is the uncertainty in the X^Zi/ cross-feed. 
It is estimated by separately varying the fraction of B —>■ ttIu and B —>■ plv decays (which 
are expected to dominate in the high region) within their respective experimental un- 



7 



TABLE II: The fitted signal yield, the selection efficiency for signal candidates, the branching 
fraction and the goodness of fit (estimated by the divided by the number of degrees of freedom) . 
For each fit (using a given form- factor model), the error quoted on the signal yield and the branching 
fraction is statistical only. For the average, the first error is statistical, and the second is the spread 
around the central value. 



form- factor model 


N{B+ ujI+u) 


ee + C/i 




xVndf 


ISGW2 


383 ± 118 


5.0% 


1.0 ±0.3 


1.05 


UKQCD [9] 


384 ± 116 


4.2% 


1.2 ±0.4 


1.08 


LCSR \W] 


473 ± 141 


3.8% 


1.7 ±0.5 


1.04 


average 


414 ± 125 ± 42 




1.3 ±0.4 ±0.3 





certainties [2J|. The relative fractions of charged and neutral modes are constrained using 
isospin symmetry. For the B — >■ plv mode, we also consider the form-factor model depen- 
dence of the cross- feed 0, 0|. To estimate the uncertainy in the cross-feed from other 
B Xulv decays, the fit is repeated modeling this component once with ISGW2, and once 
with the De Fazio- Neubert model. Half of the difference between these two cases is assigned 
as a systematic uncertainty. The next-to-largest component is the uncertainty in the neu- 
trino reconstruction, track finding and cluster finding efficiency. While the latter two are 
uncorrelated, they are treated as fully correlated with the former. The XJ,v cross-feed uncer- 
tainty is estimated by varying the fractions of -B ^ D*lu^ B Dlu and B D** / D^*^nlu 
in S — >■ XJi' within ±10%, ±10% and ±30%, respectively, and summing the individual 
variations in quadrature. 

In summary, we have measured the B^ — > ul^u branching fraction to be (1.3±0.4(stat) ± 
0.2{syst) ± 0.3{model)) x 10~^, based on 414 ± 125 signal events. This is the first evidence 
for this decay. Assuming the quark model relation F(i?'' p~l^u) = 2F(i?+ ul^u), our 
measurement agrees with measurements of the decay 5° p^l^u (sl, 0|- 
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TABLE III: Contributions to the systematic uncertainty. The size of each contribution is given as 
percentage of the branching fraction. 





B Tilv cross-feed 


2.2% 


B plv cross-feed 


14.8% 


other B — > X^lv cross-feed 


1.4% 


(sum) 


15.1% 


neutrino reconstruction 


4% 


charged track finding (Z,7r+,7r^) 


3% 


cluster finding (vr*^) 


4% 


(sum) 


9% 


Xclv cross-feed 


2.7% 


lepton identification 


3.0% 


number of BB 


0.6% 


uj TT+vr^TT^ branching fraction [24] 


0.8% 


total systematic uncertainty 


18.1% 
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